Limited associations between MHC diversity and reproductive success in a bird species with biparental care.

The selective pressure from pathogens on individuals can have direct consequences on reproduction. Genes from the major histocompatibility complex (MHC) are central to the vertebrate adaptive immune system and pathogen resistance. In species with biparental care, each sex has distinct reproductive roles and levels of investment, and due to a trade-off with immunity, one can expect different selective regimes acting upon the MHC of each parent. Here, we addressed whether couples combine each other's variation at MHC loci to increase their breeding success. Specifically, we used a 23-year dataset from a barn owl population (Tyto alba) to understand how MHC class Iα and IIß functional divergence and supertypes of each parent were associated with clutch size and fledging success. We did not detect associations between MHC diversity and supertypes with the clutch size or with the fledging success. In addition, to understand the relative contribution from the MHC of the genetic parents and the social parents, we analyzed the fledging success using only a cross-fostered dataset. We found several associations of weak-to-moderate effect sizes between the father's MHC and fledging success: (i) lower MHC-Iα divergence in the genetic father increases fledging success, which might improve paternal care during incubation, and (ii) one and two MHC-IIß DAB2 supertypes in the social father decrease and increase, respectively, fledging success, which may affect the paternal care after hatching. Furthermore, fledging success increased when both parents did not carry MHC-IIß DAB1 supertype 2, which could suggest conditional effects of this supertype. Although our study relied on a substantial dataset, we showed that the associations between MHC diversity and reproductive success remain scarce and of complex interpretation in the barn owl. Moreover, our results highlighted the need to incorporate more than one proxy of reproductive success and several MHC classes to capture more complex associations.

The genomic architecture of continuous plumage colour variation in the European barn owl (Tyto alba).

The maintenance of colour variation in wild populations has long fascinated evolutionary biologists, although most studies have focused on discrete traits exhibiting rather simple inheritance patterns and genetic architectures. However, the study of continuous colour traits and their potentially oligo- or polygenic genetic bases remains rare in wild populations. We studied the genetics of the continuously varying white-to-rufous plumage coloration of the European barn owl (Tyto alba) using a genome-wide association approach on the whole-genome data of 75 individuals. We confirmed a mutation at the melanocortin-1-receptor gene (MC1R) is involved in the coloration and identified two new regions, located in super-scaffolds 9 and 42. The combination of the three regions explains most of the colour variation (80.37%, 95% credible interval 58.45-100%). One discovered region, located in the sex chromosome, differs between the most extreme colorations in owls sharing a specific MC1R genotype. This region may play a role in the colour sex dimorphism of this species, possibly in interaction with the autosomal MC1R. We thus provide insights into the genetic architecture of continuous colour variation, pointing to an oligogenic basis with potential epistatic effects among loci that should aid future studies understanding how continuous colour variation is maintained in nature.

Plastic responses to warmer climates: a semi-natural experiment on lizard populations.

Facing warming environments, species can exhibit plastic or microevolutionary changes in their thermal physiology to adapt to novel climates. Here, using semi-natural mesocosms, we experimentally investigated over two successive years whether a 2°C-warmer climate produces selective and inter- and intragenerational plastic changes in the thermal traits (preferred temperature and dorsal coloration) of the lizard Zootoca vivipara. In a warmer climate, the dorsal darkness, dorsal contrast, and preferred temperature of adults plastically decreased and covariances between these traits were disrupted. While selection gradients were overall weak, selection gradients for darkness were slightly different between climates and in the opposite direction to plastic changes. Contrary to adults, male juveniles were darker in warmer climates either through plasticity or selection and this effect was strengthened by intergenerational plasticity when juveniles' mothers also experienced warmer climates. While the plastic changes in adult thermal traits alleviate the immediate overheating costs of warming, its opposite direction to selective gradients and to juveniles' phenotypic responses may slow down evolutionary shifts toward phenotypes that are better adapted to future climates. Our study demonstrates the importance of considering inter- and intragenerational plasticity along with selective processes to better understand adaptation and population dynamics in light of climate change.

Investigating the genetic basis of vertebrate dispersal combining RNA-seq, RAD-seq and quantitative genetics.

Although animal dispersal is known to play key roles in ecological and evolutionary processes such as colonization, population extinction and local adaptation, little is known about its genetic basis, particularly in vertebrates. Untapping the genetic basis of dispersal should deepen our understanding of how dispersal behaviour evolves, the molecular mechanisms that regulate it and link it to other phenotypic aspects in order to form the so-called dispersal syndromes. Here, we comprehensively combined quantitative genetics, genome-wide sequencing and transcriptome sequencing to investigate the genetic basis of natal dispersal in a known ecological and evolutionary model of vertebrate dispersal: the common lizard, Zootoca vivipara. Our study supports the heritability of dispersal in semi-natural populations, with less variation attributable to maternal and natal environment effects. In addition, we found an association between natal dispersal and both variation in the carbonic anhydrase (CA10) gene, and in the expression of several genes (TGFB2, SLC6A4, NOS1) involved in central nervous system functioning. These findings suggest that neurotransmitters (serotonin and nitric oxide) are involved in the regulation of dispersal and shaping dispersal syndromes. Several genes from the circadian clock (CRY2, KCTD21) were also differentially expressed between disperser and resident lizards, supporting that the circadian rhythm, known to be involved in long-distance migration in other taxa, might affect dispersal as well. Since neuronal and circadian pathways are relatively well conserved across vertebrates, our results are likely to be generalisable, and we therefore encourage future studies to further investigate the role of these pathways in shaping dispersal in vertebrates.

Molecular evolution of the proopiomelanocortin system in Barn owl species.

Examination of genetic polymorphisms in outbred wild-living species provides insights into the evolution of complex systems. In higher vertebrates, the proopiomelanocortin (POMC) precursor gives rise to α-, ß-, and γ-melanocyte-stimulating hormones (MSH), which are involved in numerous physiological aspects. Genetic defects in POMC are linked to metabolic disorders in humans and animals. In the present study, we undertook an evolutionary genetic approach complemented with biochemistry to investigate the functional consequences of genetic polymorphisms in the POMC system of free-living outbred barn owl species (family Tytonidae) at the molecular level. Our phylogenetic studies revealed a striking correlation between a loss-of-function H9P mutation in the ß-MSH receptor-binding motif and an extension of a poly-serine stretch in γ3-MSH to ≥7 residues that arose in the barn owl group 6-8 MYA ago. We found that extension of the poly-serine stretches in the γ-MSH locus affects POMC precursor processing, increasing γ3-MSH production at the expense of γ2-MSH and resulting in an overall reduction of γ-MSH signaling, which may be part of a negative feedback mechanism. Extension of the γ3-MSH poly-serine stretches ≥7 further markedly increases peptide hormone stability in plasma, which is conserved in humans, and is likely relevant to its endocrine function. In sum, our phylogenetic analysis of POMC in wild living owls uncovered a H9P ß-MSH mutation subsequent to serine extension in γ3-MSH to 7 residues, which was then followed by further serine extension. The linked MSH mutations highlight the genetic plasticity enabled by the modular design of the POMC gene.

New genome assembly of the barn owl (Tyto alba alba).

New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single-Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations.

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

Differential fitness effects of moonlight on plumage colour morphs in barn owls.

The Moon cycle exposes nocturnal life to variation in environmental light. However, whether moonlight shapes the fitness of nocturnal species with distinct colour variants remains unknown. Combining data from long-term monitoring, high-resolution global positioning system tracking and experiments using prey, we show that barn owls (Tyto alba) with distinct plumage colourations are differently affected by moonlight. The reddest owls are less successful at hunting and providing food to their offspring during moonlit nights, which associates with lower body mass and lower survival of the youngest nestlings and with female mates starting to lay eggs at low moonlight levels. Although moonlight should make white owls more conspicuous to prey, it either positively affects or does not affect the hunting and fitness of the whitest owls. We experimentally show that, under full-moon conditions, white plumage triggers longer freezing times in prey, which should facilitate prey catchability. We propose that the barn owl's white plumage, a rare trait among nocturnal predators, exploits the known aversion of rodents to bright light, explaining why, counterintuitively, moonlight has a lesser impact on the whitest owls. Our study provides evidence for the long-suspected influence of the Moon on the evolution of colouration in nocturnal species, highlighting the importance of colour in nocturnal ecosystems.

Expression of glucocorticoid and mineralocorticoid receptor genes co-varies with a stress-related colour signal in barn owls.

Glucocorticoid hormones are important intermediates between an organism and its environment. They enable an organism to adjust its behavioural and physiological processes in response to environmental changes by binding to mineralocorticoid receptors (MR) and glucocorticoid receptors (GR) expressed in many tissues, including the integument. The regulation of glucocorticoids co-varies with melanin-based colouration in numerous species, an association that might result from pleiotropic effects of genes in the melanocortin system and evolve within a signalling context. Most studies have focused on the circulating levels of glucocorticoids disregarding the receptors that mediate their action, and that might partly account for the covariation between the regulation of stress and melanin-based colouration. We investigated the association of the expression levels of GR and MR genes with melanin-based colouration in the growing feathers of nestling barn owls (Tyto alba). We also explored the association between GR and MR expression levels and the expression of genes related to the melanocortin system and melanogenesis to better understand the origin of the link between the expression of receptors to which corticosterone binds and melanin-based colouration. Nestling barn owls displaying larger eumelanic black feather spots expressed GR and MR at lower levels than smaller-spotted individuals. However, we found that the expression of the GR and MR genes was positively rather than negatively correlated with the expression of genes involved in the deposition of melanin pigments at the time we sampled the nestlings. This provides mixed evidence of the association between melanin-based traits and MR and GR gene expression. The finding that the expression of GR and MR was positively associated with the expression of the PCSK2 gene (encoding one of the protein convertase responsible for the production of hormone peptide ACTH and α-MSH) suggests that the melanocortin system may be implicated in the establishment of the covariation between melanin-based colour and the expression of receptors to which glucocorticoids bind. However, further studies investigating the expression of melanin-based traits with stress-related endpoints at different time points of feather development will be necessary to understand better the proximate mechanism linking melanin-based traits with stress.

Baseline and stress-induced corticosterone levels are heritable and genetically correlated in a barn owl population.

The hypothalamic-pituitary-adrenal (HPA) axis is responsible for the regulation of corticosterone, a hormone that is essential in the mediation of energy allocation and physiological stress. As a continuous source of challenge and stress for organisms, the environment has promoted the evolution of physiological adaptations and led to a great variation in corticosterone profiles within or among individuals, populations and species. In order to evolve via natural selection, corticosterone levels do not only depend on the strength of selection exerted on them, but also on the extent to which the regulation of corticosterone is heritable. Nevertheless, the heritability of corticosterone profiles in wild populations is still poorly understood. In this study, we estimated the heritability of baseline and stress-induced corticosterone levels in barn owl (Tyto alba) nestlings from 8 years of data, using a multivariate animal model based on a behavioural pedigree. We found that baseline and stress-induced corticosterone levels are strongly genetically correlated (r = 0.68-0.80) and that the heritability of stress-induced corticosterone levels (h2 = 0.24-0.33) was moderate and similar to the heritability of baseline corticosterone levels (h2 = 0.19-0.30). These findings suggest that the regulation of stress-induced corticosterone and baseline levels evolves at a similar pace when selection acts with the same intensity on both traits and that contrary to previous studies, the evolution of baseline and stress-induced level is interdependent in barn owls, as they may be strongly genetically correlated.

Toward Understanding the Repeated Occurrence of Associations between Melanin-Based Coloration and Multiple Phenotypes.

Melanin is the most widespread pigment in organisms. Melanin-based coloration has been repeatedly observed to be associated with the same traits and in the same direction in different vertebrate and insect species. However, whether any factors that are common to different taxa account for the repeated evolution of melanin-phenotype associations remains unclear. We propose to approach this question from the perspective of convergent and parallel evolution to clarify to what extent different species have evolved the same associations owing to a shared genetic basis and being subjected to similar selective pressures. Our current understanding of the genetic basis of melanin-phenotype associations allows for both convergent and parallel evolution, but this understanding is still limited. Further research is needed to clarify the generality and interdependencies of the different proposed mechanisms (supergenes, pleiotropy based on hormones, or neural crest cells). The general ecological scenarios whereby melanin-based coloration is under selection-protection from ultraviolet radiation, thermoregulation in cold environments, or as a signal of social status-offer a good opportunity to study how melanin-phenotype associations evolve. Reviewing these scenarios shows that some traits associated with melanin-based coloration might be selected together with coloration by also favoring adaptation but that other associated traits might impede adaptation, which may be indicative of genetic constraints. We therefore encourage further research on the relative roles that selection and genetic constraints play in shaping multiple melanin-phenotype associations. Placed into a phylogenetic context, this will help clarify to what extent these associations result from convergent or parallel evolutionary processes and why melanin-phenotype associations are so common across the tree of life.

Selection on the Major Color Gene Melanocortin-1-Receptor Shaped the Evolution of the Melanocortin System Genes.

Modular genetic systems and networks have complex evolutionary histories shaped by selection acting on single genes as well as on their integrated function within the network. However, uncovering molecular coevolution requires the detection of coevolving sites in sequences. Detailed knowledge of the functions of each gene in the system is also necessary to identify the selective agents driving coevolution. Using recently developed computational tools, we investigated the effect of positive selection on the coevolution of ten major genes in the melanocortin system, responsible for multiple physiological functions and human diseases. Substitutions driven by positive selection at the melanocortin-1-receptor (MC1R) induced more coevolutionary changes on the system than positive selection on other genes in the system. Contrarily, selection on the highly pleiotropic POMC gene, which orchestrates the activation of the different melanocortin receptors, had the lowest coevolutionary influence. MC1R and possibly its main function, melanin pigmentation, seems to have influenced the evolution of the melanocortin system more than functions regulated by MC2-5Rs such as energy homeostasis, glucocorticoid-dependent stress and anti-inflammatory responses. Although replication in other regulatory systems is needed, this suggests that single functional aspects of a genetic network or system can be of higher importance than others in shaping coevolution among the genes that integrate it.

Beyond mean allelic effects: A locus at the major color gene MC1R associates also with differing levels of phenotypic and genetic (co)variance for coloration in barn owls.

The mean phenotypic effects of a discovered variant help to predict major aspects of the evolution and inheritance of a phenotype. However, differences in the phenotypic variance associated to distinct genotypes are often overlooked despite being suggestive of processes that largely influence phenotypic evolution, such as interactions between the genotypes with the environment or the genetic background. We present empirical evidence for a mutation at the melanocortin-1-receptor gene, a major vertebrate coloration gene, affecting phenotypic variance in the barn owl, Tyto alba. The white MC1R allele, which associates with whiter plumage coloration, also associates with a pronounced phenotypic and additive genetic variance for distinct color traits. Contrarily, the rufous allele, associated with a rufous coloration, relates to a lower phenotypic and additive genetic variance, suggesting that this allele may be epistatic over other color loci. Variance differences between genotypes entailed differences in the strength of phenotypic and genetic associations between color traits, suggesting that differences in variance also alter the level of integration between traits. This study highlights that addressing variance differences of genotypes in wild populations provides interesting new insights into the evolutionary mechanisms and the genetic architecture underlying the phenotype.

Genomics of coloration in natural animal populations.

Animal coloration has traditionally been the target of genetic and evolutionary studies. However, until very recently, the study of the genetic basis of animal coloration has been mainly restricted to model species, whereas research on non-model species has been either neglected or mainly based on candidate approaches, and thereby limited by the knowledge obtained in model species. Recent high-throughput sequencing technologies allow us to overcome previous limitations, and open new avenues to study the genetic basis of animal coloration in a broader number of species and colour traits, and to address the general relevance of different genetic structures and their implications for the evolution of colour. In this review, we highlight aspects where genome-wide studies could be of major utility to fill in the gaps in our understanding of the biology and evolution of animal coloration. The new genomic approaches have been promptly adopted to study animal coloration although substantial work is still needed to consider a larger range of species and colour traits, such as those exhibiting continuous variation or based on reflective structures. We argue that a robust advancement in the study of animal coloration will also require large efforts to validate the functional role of the genes and variants discovered using genome-wide tools.This article is part of the themed issue 'Animal coloration: production, perception, function and application'.

MC1R variants affect the expression of melanocortin and melanogenic genes and the association between melanocortin genes and coloration.

The melanocortin-1 receptor (MC1R) gene influences coloration by altering the expression of genes acting downstream in the melanin synthesis. MC1R belongs to the melanocortin system, a genetic network coding for the ligands that regulate MC1R and other melanocortin receptors controlling different physiological and behavioural traits. The impact of MC1R variants on these regulatory melanocortin genes was never considered, even though MC1R mutations could alter the influence of these genes on coloration (e.g. by decreasing MC1R response to melanocortin ligands). Using barn owl growing feathers, we investigated the differences between MC1R genotypes in the (co)expression of six melanocortin and nine melanogenic-related genes and in the association between melanocortin gene expression and phenotype (feather pheomelanin content). Compared to the MC1R rufous allele, responsible for reddish coloration, the white allele was not only associated with an expected lower expression of melanogenic-related genes (TYR, TYRP1, OCA2, SLC45A2, KIT, DCT) but also with a lower MC1R expression and a higher expression of ASIP, the MC1R antagonist. More importantly, the expression of PCSK2, responsible for the maturation of the MC1R agonist, α-melanocyte-stimulating hormone, was positively related to pheomelanin content in MC1R white homozygotes but not in individuals carrying the MC1R rufous allele. These findings indicate that MC1R mutations not only alter the expression of melanogenic-related genes but also the association between coloration and the expression of melanocortin genes upstream of MC1R. This suggests that MC1R mutations can modulate the regulation of coloration by the pleiotropic melanocortin genes, potentially decoupling the often-observed associations between coloration and other phenotypes.

Inter-class competition in stage-structured populations: effects of adult density on life-history traits of adult and juvenile common lizards.

Ecological and evolutionary processes in natural populations are largely influenced by the population's stage-structure. Commonly, different classes have different competitive abilities, e.g., due to differences in body size, suggesting that inter-class competition may be important and largely asymmetric. However, experimental evidence states that inter-class competition, which is important, is rare and restricted to marine fish. Here, we manipulated the adult density in six semi-natural populations of the European common lizard, Zootoca vivipara, while holding juvenile density constant. Adult density affected juveniles, but not adults, in line with inter-class competition. High adult density led to lower juvenile survival and growth before hibernation. In contrast, juvenile survival after hibernation was higher in populations with high adult density, pointing to relaxed inter-class competition. As a result, annual survival was not affected by adult density, showing that differences in pre- and post-hibernation survival balanced each other out. The intensity of inter-class competition affected reproduction, performance, and body size in juveniles. Path analyses unravelled direct treatment effects on early growth (pre-hibernation) and no direct treatment effects on the parameters measured after hibernation. This points to allometry of treatment-induced differences in early growth, and it suggests that inter-class competition mainly affects the early growth of the competitively inferior class and thereby their future performance and reproduction. These results are in contrast with previous findings and, together with results in marine fish, suggest that the strength and direction of density dependence may depend on the degree of inter-class competition, and thus on the availability of resources used by the competing classes.

Effect of the MC1R gene on sexual dimorphism in melanin-based colorations.

Variants of the melanocortin-1 receptor (MC1R) gene result in abrupt, naturally selected colour morphs. These genetic variants may differentially affect sexual dimorphism if one morph is naturally selected in the two sexes but another morph is naturally or sexually selected only in one of the two sexes (e.g. to confer camouflage in reproductive females or confer mating advantage in males). Therefore, the balance between natural and sexual selections can differ between MC1R variants, as suggest studies showing interspecific correlations between sexual dimorphism and the rate of nonsynonymous vs. synonymous amino acid substitutions at the MC1R. Surprisingly, how MC1R is related to within-species sexual dimorphism, and thereby to sex-specific selection, has not yet been investigated. We tackled this issue in the barn owl (Tyto alba), a species showing pronounced variation in the degree of reddish pheomelanin-based coloration and in the number and size of black feather spots. We found that a valine (V)-to-isoleucine (I) substitution at position 126 explains up to 30% of the variation in the three melanin-based colour traits and in feather melanin content. Interestingly, MC1R genotypes also differed in the degree of sexual colour dimorphism, with individuals homozygous for the II MC1R variant being 2 times redder and 2.5 times less sexually dimorphic than homozygous individuals for the VV MC1R variant. These findings support that MC1R interacts with the expression of sexual dimorphism and suggest that a gene with major phenotypic effects and weakly influenced by variation in body condition can participate in sex-specific selection processes.

Cumulative frequency-dependent selective episodes allow for rapid morph cycles and rock-paper-scissors dynamics in species with overlapping generations.

Rock-paper-scissors (RPS) dynamics, which maintain genetic polymorphisms over time through negative frequency-dependent (FD) selection, can evolve in short-lived species with no generational overlap, where they produce rapid morph frequency cycles. However, most species have overlapping generations and thus, rapid RPS dynamics are thought to require stronger FD selection, the existence of which yet needs to be proved. Here, we experimentally demonstrate that two cumulative selective episodes, FD sexual selection reinforced by FD selection on offspring survival, generate sufficiently strong selection to generate rapid morph frequency cycles in the European common lizard Zootoca vivipara, a multi-annual species with major generational overlap. These findings show that the conditions required for the evolution of RPS games are fulfilled by almost all species exhibiting genetic polymorphisms and suggest that RPS games may be responsible for the maintenance of genetic diversity in a wide range of species.

Iridophores and not carotenoids account for chromatic variation of carotenoid-based coloration in common lizards (Lacerta vivipara).

Abstract Carotenoids typically need reflective background components to shine. Such components, iridophores, leucophores, and keratin- and collagen-derived structures, are generally assumed to show no or little environmental variability. Here, we investigate the origin of environmentally induced variation in the carotenoid-based ventral coloration of male common lizards (Lacerta vivipara) by investigating the effects of dietary carotenoids and corticosterone on both carotenoid- and background-related reflectance. We observed a general negative chromatic change that was prevented by ß-carotene supplementation. However, chromatic changes did not result from changes in carotenoid-related reflectance or skin carotenoid content but from changes in background-related reflectance that may have been mediated by vitamin A1. An in vitro experiment showed that the encountered chromatic changes most likely resulted from changes in iridophore reflectance. Our findings demonstrate that chromatic variation in carotenoid-based ornaments may not exclusively reflect differences in integumentary carotenoid content and, hence, in qualities linked to carotenoid deposition (e.g., foraging ability, immune response, or antioxidant capacity). Moreover, skin carotenoid content and carotenoid-related reflectance were related to male color polymorphism, suggesting that carotenoid-based coloration of male common lizards is a multicomponent signal, with iridophores reflecting environmental conditions and carotenoids reflecting genetically based color morphs.